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ABSTRACT 

We present spatially resolved spectroscopy of the supernova remnant (SNR) G292.0+1.8 with the Chandra 
X-ray observatory. This composite-type SNR contains the 135 ms pulsar, Jl 124-5916, recently discovered by 
Camilo and coworkers. We apply non-equilibrium ionization (NEI) models to fit the ejecta-dominated regions 
and to identify the blast wave of the supernova explosion. By comparing the derived abundances with those 
predicted from nucleosynthesis models, we estimate a progenitor mass of 3CMK) M . We derive the intrinsic 
parameters of the supernova explosion such as its energy, the age of the SNR, the blast wave velocity, and the 
swept-up mass. In the Sedov interpretation, our estimated SNR age of 2,600^00 years is close to the pulsar's 
characteristic age of 2,900 years. This confirms the pulsar/SNR association and relaxes the need for the pulsar to 
have a non-canonical value for the braking index, a large period at birth or a large transverse velocity. We discuss 
the properties of the pulsar wind nebula (PWN) in the light of the Kennel and Coroniti model and estimate the 
pulsar wind magnetization parameter. We also report the first evidence for steepening of the power law spectral 
index with increasing radius from the pulsar, a result that is expected from synchrotron losses and is reminiscent 
of Crab-like PWNe. 

Subject headings: ISM: individual (G292.0+1.8, MSH 1 1-54, PSR Jl 124-5916) — ISM: abundances — 
supernova remnants — X-rays: ISM 



1. INTRODUCTION 

Chandra has been advancing our understanding of core- 
collapse supernovae by revealing their collapsed cores and their 
outflows, the distribution of supernova ejecta, and their interac- 
tion with the interstellar medium (ISM). Increasing the sample 
of SNRs with a hybrid morphology (i.e. Composite-type SNRs) 
helps address questions related to the birth properties and evo- 
lution of neutron stars and SNRs in the ISM. 

G292.0+1.8 (MSH 11-54) offers a unique laboratory to ad- 
dress these questions because of its hybrid morphology, its 
youth, and its size. This SNR is one of three Oxygen-rich 
remnants in our Galaxy. It was first discovered in the radio 
band (Milne 1969; Shaver & Goss 1970), where it exhibits 
a centrally-filled morphology indicative of the presence of a 
PWN. Einstein and EXOSAT X-ray observations classified the 
remnant as a type-II explosion of a massive star (Hughes & 
Singh 1994). Chandra observations led to the discovery of the 
long-sought PWN (Hughes et al. 2001; Safi-Harb & Gonzalez 

2001) surrounding a candidate pulsar. An imaging study of the 
SNR with Chandra (Safi-Harb & Gonzalez 2001; Park et al. 

2002) traced the distribution of the ejecta. The remnant has a 
morphology similar to Cas A on small scales (Hwang, Holt, & 
Petre 2001), however, no Fe-K line-emission was detected. 

The discovery of the PWN in G292.0+1.8 prompted a deep 
search for radio pulsations and led to the discovery of a 135 ms 
pulsar coincident with the Chandra X-ray point source (Camilo 
et al. 2002a). The pulsar has a derived spin-down luminosity 
E of 1.2 x 10 37 erg s , a surface magnetic dipole field strength 
of 1 x 10 13 G, and a characteristic age of 2,900 yrs. To date, the 
only estimate of the age of this SNR is from optical studies of 
the high-velocity material coincident with the PWN (Murdin & 
Clark 1979). The age of 1,600 years is a factor of 1.8 less than 
the characteristic age of the pulsar. 

In this Letter, we present the first detailed spatially resolved 
spectroscopy of the Chandra observation in order to infer the 
intrinsic parameters of the supernova explosion and its associ- 



ated pulsar and PWN. In §2, we briefly summarize the obser- 
vation and present our sub-arcminute spatially resolved spec- 
troscopy of the SNR and the embedded PWN. In §3, we discuss 
our results and infer the energetics of the supernova explosion, 
the mass of the progenitor, the SNR's age, the pulsar's trans- 
verse velocity and the pulsar wind parameter. We also report 
the first evidence for spectral softening of the power law index 
of the PWN with increasing distance from the pulsar. 

2. SUB-ARCMINUTE SPATIALLY RESOLVED SPECTROSCOPY 

G292.0+1.8 was observed with the S3-chip of the Advanced 
CCD Imaging Spectrometer (ACIS) on board Chandra on 
March 11, 2000, at a CCD temperature of -120"C. We cor- 
rected for Charge Transfer Inefficiency (CTI) by applying the 
cticorrectit tool to the event 1 raw data (Townsley et al. 2000). 
A new event 2 file was then obtained by screening the data 
using standard Ciao 2.2 routines. The resulting exposure time 
was 38 ksec. 

2.1. SNR 

To perform spatially resolved spectroscopy of the SNR, we 
selected the 23 bright regions shown in Fig. 1 (left). In this 
image, the broadband emission from the SNR was divided into 
soft (0.3-1.15 keV, red), medium (1.15-2.15 keV, green) and 
hard (2.15-10.0 keV, blue) contributions, each smoothed with 
a gaussian with ct=1."5 . The PWN is clearly visible in blue 
near the center of the SNR. We fitted the spectra of the rem- 
nant in the 0.5-8.0 keV range using the non-equilibrium ion- 
ization (NEI) model, vpshock (Borkowski, Lyerly, & Reynolds 
2001) in version 11.0.1 of the XSPEC package (Arnaud 1996). 
This model is appropriate for modeling plane-parallel shocks in 
young SNRs whose plasma has not reached ionization equilib- 
rium. It is characterized by the shock temperature, kT, and the 
ionization timescale, n e t, of the plasma. Here, n e is the post- 
shock electron density and t is the time since the passage of the 
shock, t can then be used as a lower estimate on the SNR age. 
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Since the remnant is ejecta-dominated, we allowed the abun- 
dances to vary. 

One-component NEI model fits yielded an average temper- 
ature for the selected regions of 0.9±0.3 keV, with minimum 
and maximum values of 0.6 and 1 .7 keV, respectively. The ion- 
ization timescale varied across the SNR from ^10 n to 10 13 
cm" 3 s. These values indicate that the emission from the ejecta 
contains a wide range of temperatures and ionization states. 
The highest ionization state values are associated with regions 
1, 3, and 9 (see Fig. 1, left), indicating that the plasma has 
reached ionization equilibrium in the southeast. Relative abun- 
dances were well constrained in these models and were used to 
derive the abundance ratios discussed in §3.2. 

We subsequently attempted two-component model fits. This 
was motivated by the expectation of a high-temperature and 
a low-temperature plasma associated with the supernova blast 
wave and reverse shock, respectively (see e.g. Safi-Harb et al. 
2000). While the overall fits were slightly improved, the abun- 
dances could not be well constrained due to the large number 
of free parameters. The derived average temperatures for the 
hot and cool components were 1.05±0.34 keV and 0.37±0.18 
keV, respectively. These are lower than 1.64^2 keV, as re- 
ported by Hughes & Singh (1994) from their one-component 
NEI model fit to the Einstein and EXOSAT data. The discrep- 
ancy is probably due to the contamination by the previously 
unresolved PWN. 

2.2. PWN 

Chandra observations of PWNe have shown that 1) PWNe 
harbor torus-like and jet-like structures indicating equatorial 
and polar outflows, and 2) the power law photon index steep- 
ens away from the pulsar, as expected from synchrotron and 
expansion losses (see e.g. Safi-Harb 2002 for a review). By 
examining this PWN, we find evidence of an arcsecond-scale 
east- west elongation and an extension to the south (see Fig 1, 
right), hinting at structures associated with the deposition of the 
pulsar's wind energy into its surroundings. To search for spec- 
tral variations, we divided the inner 45" of the PWN into four 
concentric rings to obtain enough counts to extract a spectrum. 
As reported by Hughes et al. (2001), the spectrum of the PWN 
is heavily contaminated by thermal emission from the SNR. In 
order to minimize this contamination, we used a nearby region 
as background and restricted the fit to the 1.5-7.0 keV range, 
thus reducing the contamination by Oxygen, Neon, and Mag- 
nesium ejecta. Using a power-law model, we find that the pho- 
ton index (T) steepens away from the pulsar (see Table 1). This 
result is expected from synchrotron losses and is reminiscent of 
Crab-like nebulae (e.g. Safi-Harb et al. 2001). 

3. DISCUSSION 
3.1. Distance 

From our one- and two-component fits, we derive an aver- 
age value of Nh=(0.5 ± 0.1) x 10 22 cm" 2 for the remnant. We 
can use this value to verify previous estimates for the distance 
to G292.0+1.8. The extinction per unit distance in the direc- 
tion of the remnant can be estimated from the contour diagrams 
given by Lucke (1978): E B - V /D~0.2 mag kpc" 1 . Using the 
relation 7Y///Av=l. 79 x 10 21 cm" 2 mag" 1 , which translates into 
7Y H /£s-v=5.55xl0 21 cm" 2 mag" 1 (Predehl & Schmitt 1995), 
we derive a distance of 3.6-5.5 kpc. A distance of 4.8 kpc, as 
used in the previous studies, is then a good approximation and 
is adopted in the following discussion. 



3.2. Supernova Parameters 

We subsequently use our NEI model fits to derive the su- 
pernova parameters. In the standard picture of X-ray emis- 
sion from young SNRs, the soft component arises from shocked 
ejecta and the hard component is usually attributed to the blast 
wave shocking the surrounding medium. Our spatially resolved 
spectroscopy indicates that most of the regions in the SNR 
(except for the PWN, the outermost regions and the eastern 
'belt' -like filament near the SNR center) are well fitted with 
a one-component NEI (vpshock) model with high metal abun- 
dances (with respect to solar), indicating that these knots are the 
shocked ejecta associated with the explosion of a massive star. 
Nucleosynthesis yield as a function of the progenitor's mass has 
been numerically computed by Woosley and Weaver (1995) and 
several others. In order to derive the progenitor's mass, we use 
the abundance ratios from our NEI model fits in the 0.5-8.0 
keV band to the ejecta-dominated regions. In Fig. 2, we plot 
the average values for these ratios along with their rms scat- 
ter. Over-plotted are the predicted ratios from nucleosynthesis 
models for progenitors of 25, 30, 35 and 40 M Q . We then con- 
clude that the measured abundance ratios are consistent with a 
progenitor's mass of 30^-0 M©. 

Regions 20 to 23 selected in our analysis are of special in- 
terest (see Fig. 1, left). Region 20 is part of a set of sharp 
circumferential filaments surrounding the bulk of the ejecta- 
dominated X-ray emission. Regions 21, 22 and 23 are part of 
an outer faint ridge of emission delimiting the SNR. We be- 
lieve that these latter regions indicate the location of the blast 
wave shocking the surrounding medium for the following rea- 
sons: 1) as shown in Table 2, the abundances derived for re- 
gion 21 using one-component models are consistent with ISM 
abundances (sub-solar) and are the lowest ones throughout the 
remnant, while those derived for region 20 are contaminated 
by ejecta; 2) the outermost emission is absent in the equiva- 
lent width maps which reveal the ejecta (Park et al. 2002); and 
3) even two-component model fits to regions 2 1-23 require sub- 
solar abundances, which is not the case for region 20, where 
both components require above-solar abundances. 

In order to estimate the SNR age, we will adopt our more re- 
liable one-component fits to regions 21-23 in the 0.5-8.0 keV 
band. The corresponding SNR size would be r^4.2'=5.9 pc 
(at 4.8 kpc). In the following example, we will use the pa- 
rameters derived for region 21 shown in Table 2. The ob- 
served emission measure (EM) of 0.09±0.02 (after scaling it to 
the SNR volume) corresponds to: / n e n H dV=10 14 (47rD 2 )(EM) 
= (2.5±0.6)xl0 58 cm" 3 . If the emission volume V=fV tot , 
where / is the filling factor, we find an upstream density 
no = (0.23±0.03) f~ l l 2 cm" 3 . From our fit parameter not, 
obtained from the vpshock model, we find a shock age of 
f=(2,500^7o[J°)/ 1 ' /2 years. Similarly, we derive the parameters 
for regions 22 and 23 and summarize them in Table 3. For a rea- 
sonable filling factor of 0.25, the shock age varies from 1,250- 
4,450 years. The discrepancy between these ages suggests that 
the northeastern regions are shocked at earlier times that the 
western regions, probably due to an inhomogeneous medium. 
This complication makes the determination of the SNR param- 
eters in this manner uncertain. 

In order to estimate the global properties of the SNR inde- 
pendently of the small-scale density variations, we attempted 
to fit the spectrum of the whole remnant. We found this 
to be extremely difficult using the entire 0.5-8.0 keV band 
and various combination models with up to three components. 
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This difficulty is due to the large number of NEI components 
needed to account for the variations in temperature and ion- 
ization time-scale of the ejecta, as indicated in our spatially 
resolved spectroscopy described above. However, restricting 
the fit to the 2.0-8.0 keV range (to remove the contamination 
from O, Ne, Mg and Si ejecta) and excluding the pulsar and 
its PWN, we found that a two-component model (vpshock + 
sedov) provides an acceptable fit for the overall spectrum with 
X^(dof)=1.17(197); sedov is an NEI model that accounts for the 
hot thermal emission from the SNR interior and is based on the 
Sedov dynamics (Borkowski et al. 2001). The vpshock model 
was used to characterize the emission from the Sulfur ejecta, 
while the sedov model was used to characterize the blast wave 
component (abundances of 2/3 solar with equal electron and 
ion temperatures). In Table 4, we summarize the SNR param- 
eters estimated from using equations (4a-4f) in Hamilton et al. 
(1983). The derived tabulated value of «o is consistent within 
the various uncertainties with that derived from the observed 
emission measure. The SNR age, f(Sedov)=2,400-2,850 yrs, 
is comparable to the average value derived from regions 21-23 
above. This age represents a more reliable estimate for the age 
of the remnant since the contamination from ejecta has been 
minimized and the effect of an inhomogeneous medium has 
been averaged. In addition, this age is in close agreement with 
the pulsar's characteristic age of 2,900 years. This agreement 
strengthens the pulsar/SNR association and relaxes the need for 
invoking a very slow initial spin period or a non-canonical value 
for the braking index of the pulsar (Camilo et al. 2002a). 

3.3. Pulsar Wind Nebula Parameters 

As noted by Hughes et al. (2001), the pulsar and its PWN 
are not located at the geometrical center of the remnant. The 
offset was attributed to either a high transverse velocity for the 
pulsar or slower expansion of the SNR toward the south-east. 
From the offset (~1') they calculate a transverse velocity for 
the pulsar of ^770(ii4.8/fi6oo) km s" 1 , where dn.% is the dis- 
tance in units of 4.8 kpc and fieoo is the age of the remnant 
in units of 1,600 yrs. Using our new estimate of 2,600^00 Y rs 
for the age we obtain a transverse velocity for the pulsar of 
470±40 c/4.8 km s" 1 . This value is in excellent agreement with 
recent estimates of the mean birth velocity of pulsars (450±90 
km s" 1 , Lyne & Lorimer 1994). The nearly east-west exten- 
sion and jet-like structure extending to the south (see Fig. 1, 
right) could be indicative of a toroidal and collimated outflows 
reminiscent of Crab-like PWNe. 

In the Kennel & Coroniti (1984) model, the synchrotron neb- 
ula is produced by the confinement of a relativistic wind in- 
jected by the pulsar into its surroundings. The radius of the 
standing shock, R s , is estimated by equating the pressure of 
the pulsar's wind, EI{AircR 2 s ), with the pressure in the nebula, 
P„=B 2 1 IAn. B n , the nebular magnetic field, can be approximated 
from its equipartition value, B n , eq . Using the average lumi- 
nosity and spectral index for the PWN, we estimate a value 
of B n . eq -1 x 10" 4 G. With the measured E=1.2x 10 37 ergs s _1 
(Camilo et al. 2002a), we obtain a shock radius of /? s ^3" at 4.8 
kpc. Alternatively, we can estimate the shock radius using the 
thermal pressure in the SNR interior as a confining pressure. 
This average pressure can be estimated as 2n e kT <~ 6.9 xlO" 9 
erg cm" 3 ; which yields a shock radius of 1" . The inner cir- 
cle shown in Fig. 1 (right) shows the estimated upper value of 
R s . An interesting parameter than can also be derived from 
Chandra studies of PWNe is the pulsar wind magnetization 
parameter, a, defined as the ratio of the Poynting flux to the 
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particles flux, which can be estimated as a ~ (R s /R„) 2 . Using 
a shock radius /?,~l"-3" and a nebular radius /?„^30"-45", 
we find that cr~5 x lO^-O.Ol, indicating a particle dominated 
wind. 

The properties of the PWN can thus be derived from the Ken- 
nel & Coroniti model. The value derived for the pulsar wind 
magnetization parameter in this manner is comparable to that 
derived for other PWNe of similar properties (Safi-Harb 2002). 
It is worth noting that the properties of the pulsar and the PWN 
in G292.0+1.8 are strikingly similar to those of the pulsar dis- 
covered in the Crab's cousin, G54. 1+0.2 (Camilo et al. 2002b, 
Lu et al. 2002). The lack of ejecta and a shell in the latter SNR 
(and in a handful other plerions of comparable age) indicates 
that the classification of plerionic-type SNRs is to a large ex- 
tent linked to the nature of the progenitor star, the explosion 
energy, and the environment. 
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Table 1 

Steepening of the power law photon index of the PWN in G292.0+ 1 .8" 



Radius 


T (± 90%) 


X 2 (dof) 


< 2" (pulsar) 


1.74±0.10 


142 (140) 


2" 7" 


2.0±0.10 


122 (94) 


7" - 15" 


2.06±0.11 


129 (125) 


15" - 30" 


2.37±0.10 


202 (219) 


30" - 45" 


3.00±0.20 


204 (212) 



a Nh fixed to the best fit value for the pul- 
sar of 0.37 x 10 22 cm" 2 



Table 2 

One- and two-component NEI fits to regions 20 and 21 



Region 20 Region 21 



N H (10 21 cirT 2 ) 


4.7 ± 0.2 


4.9 ±0.3 


kT (keV) 


66 +001 


0.87±0.07 


n e t (cm -3 s) 


(7.5^) xlO 11 


(8.7±2.4)xl0 10 


O(M ) 


1.1 ± 0.1 


n OO+0.04 
U -> 8 -0.06 


Ne (M ) 


0.90 ± 0.06 


53+ 05 

u - ,J -0.06 


Mg (M ) 


63 +005 

u - o; '-0.06 


0.23 ± 0.06 


Si (M ) 


0.47 ± 0.06 


n 97+0.09 
u - z/ -0.08 


S(M ) 


1.2 ±0.3 


2+ 04 


Fe (M ) 


0.24 ± 0.02 


0.16 ±0.03 


norm 


(3.3 ±0.2)xl0" 3 


(1.4±0.2)xl0" 3 


X* (dof) 


1.9(119) 


1.2 (98) 


Afo (10 zi cm" 2 ) 


6.1 ±0.03 


6 1 +,U)2 


kT hot (keV) 


94+0 05 


95+0- 20 


n e f/, , (cm -3 s) 


(3 ±0.5)xl0 n 


(6!i't)xl0 9 


norm/,,,, 


(1.05 ±0.05) xlO" 3 


(9.1 ± 1.4) xlO" 4 


fer coo , (keV) 


0.27 ± 0.02 


0.84 ± 0.06 


« c f coo ; (cm -3 s) 


(2.7! 23 )xlO n 


(1.4 ±0.3) xlO 11 


norm coo/ 


(5.1 ±0.4)xl0" 3 


(2.8^)xl0" 4 


Xl (dof) 


1.6(116) 


1.0 (95) 



Table 3 

Derived SNR parameters from Regions 21, 22 and 23 

Region 21 Region 22 Region 23 

Ionization timescale, n e t (cm -3 s) 
Jn e n H dV (cm -3 ) 
Ambient density, no (/~ 1//2 cm" 3 ) 
Swept-up mass, M s (/~ 1/2 M ) 
Lower age limit, t (f 1 ' 2 years) 



(8.7±2.4)xl0 10 (2.1±0.6)xlO n (3.1^)xl0 u 

(2.5±0.6)xl0 58 (3.9±0.9)xl0 58 (2.5 ±0.8)x 10 58 

0.23±0.03 0.29±0.4 0.23±0.04 

6.7±0.8 8.3±0.9 6.7±1.1 

2500+ i,ioo 4 900 +2 ' 100 8 900 +5 ™° 
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Table 4 

Explosion parameters from sedov fit to hard component of the overall spectrum 



Parameter 




Value (± la) 


Emission Measure, EM 




a i 66+0.006 


Shock temperature, kT s (keV) 




0.95 ±0.15 


Ionization timescale, n e t (cm" 




(2.0!°^) xlO 11 


Shock velocity, v s (km s -1 ) 




880 ± 70 


Age, t (yrs) 




2 fSOO +25 ° 


Ambient density, no (cm -3 ) 




51 +l)15 


Explosion energy, Eq (10 51 er| 


JS) 




Swept-up mass, M s (M Q ) 




15.6^:1 




FIG. 1 . — Left: combined broadband image of G292.0+1.8; individual images were used from the soft (0.3-1.15 keV, red), medium (1.15-2.15 keV, green) and 
hard (2.15-10.0 keV, blue) bands each smoothed with a gaussian with <r=l"5. Superimposed are the regions used for our spatially resolved spectroscopy. Right: 
2.6-10.0 keV false color image of the PWN surrounding PSR Jl 124-5916. The image was smoothed using a gaussian with a = Of'25. The inner circle is at a radius 
of 3" (fi.0 and the outer one is at 45" (R n ). 
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FIG. 2. — Average abundance ratios (filled squares) and their rms scatter measured throughout the remnant. Predicted ratios from nucleosynthesis models are also 
shown. 



